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ABSTRACT 



Context. Models of stellar structure and evolution can be constrained by measuring accurate parameters of detached eclipsing binaries 
in open clusters. Multiple binary stars provide the means to determine helium abundances in these old stellar systems, and in turn, to 
improve estimates of their age. 

Aims. In the first paper of this series, we demonstrated how measurements of multiple eclipsing binaries in the old open cluster 
NGC 6791 sets tighter constraints on the properties of stellar models than has previously been possible, thereby potentially improving 
both the accuracy and precision of the cluster age. Here we add additional constraints and perform an extensive model comparison to 
determine the best estimates of the cluster age and helium content, employing as many observational constraints as possible. 
Methods. We improve our photometry and correct empirically for differential reddening effects. We then perform an extensive com- 
parison of the new colour-magnitude diagrams (CMDs) and eclipsing binary measurements to Victoria and DSEP isochrones in order 
to estimate cluster parameters. We also reanalyse a spectrum of the star 2-17 to improve [Fe/H] constraints. 

Results. We find a best estimate of the age of ~ 8.3 Gyr for NGC 6791 while demonstrating that remaining age uncertainty is 
dominated by uncertainties in the CNO abundances. The helium mass fraction is well constrained at Y = 0.30 ± 0.01 resulting in 
AK/AZ ~ 1.4 assuming that such a relation exists. During the analysis we firmly identify blue straggler stars, including the star 2-17, 
and find indications for the presence of their evolved counterparts. Our analysis supports the RGB mass-loss found from asteroseis- 
mology and we determine precisely the absolute mass of stars on the lower RGB, Mrgb = 1.15 ± 0.02Mq. This will be an important 
consistency check for the detailed asteroseismology of cluster stars. 

Conclusions. Using multiple, detached eclipsing binaries for determining stellar cluster ages, it is now possible to constrain param- 
eters of stellar models, notably the helium content, which were previously out of reach. By observing a suitable number of detached 
eclipsing binaries in several open clusters, it will be possible to calibrate the age-scale and the helium enrichment parameter A Y/AZ, 
and provide firm constraints that stellar models must reproduce. 

Key words. Open clusters: individual NGC 6791 - Stars: evolution - Stars: binaries: eclipsing - Stars: abundances - Techniques: 
spectroscopic - Techniques: photometric 



1. Introduction 

The open cluster NGC 6791 is one of the oldest and, at 
the same time, the most metal-rich open clusters known 
(Origl ia et al. 2006; ICarretta et al.1l2007l lAnthonv-Twarog et alj 
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* Based on observations carried out at the Nordic Optical Telescope 
at La Palma and ESO's VLT/UVES ESO, Paranal, Chile (75.D-0206A, 
77.D-0827A, 81.D-0091). 

** Tables 5 and 6 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



12007 : Brogaa rd et all 1201 It hereafter paper I), and is there- 
fore very important for understanding chemical evolution and 
metal-rich stars. In addition, it is a well-populated open clus- 
ter, with a (statistically) significant number of stars in all 
stages of evolution from the main sequence to the white dwarfs 
( Kin g et alJl2005l: iBedin et alJl200ll2008ah. as well as with nu- 
merous variable stars (Bruntt et al.l l2003l iMochejska et al.f 2002: 
Ide MarchietaTll2007l) . 

Despite a number of studies that have presented well- 
calibrated photo metry and high-precision colour-mag nitude dia- 
grams (CMDs) dKing et alJl2005l:IStetson et alJl2 003). the age of 
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NGC 6791 has remained very uncertain until recently because of 
correlated uncertainties in distance, reddening, and metallicity. 

It is widely appreciated that detached eclipsing binaries offer 
the possibility of determining accurate (and precise) masses and 
radii for the s ystem components, nearly indepen dent of model 
assumptions dAndersenlll99lt iTorres et aT1l2010l) . If the binary 
resides in a star cluster, and one or both of its components are 
close to the turn-off, it is possible to put tight constraints on 
the age of the system by comparing the position of the primary 
and secondary in a mass-radius (MR) diagram with theoretical 
isochrones. For stellar clusters, such an analysis has some sig- 
nificant advantages: the determination of the masses and radii 
is independent of the usual uncertainties such as reddening and 
distance. Furthermore, since the comparison to models is car- 
ried out in the MR diagram, one avoids the difficult process of 
transforming the effective temperatures and luminosities of the 
models to observed colours and magnitudes. Thus, determining 
cluster ages in the MR diagram allows a direct comparison be- 
tween observations and th eory. 

iGrundahl et al.l (120081) showed that, by applying this method 
on the eclipsing binary V20, a cluster member, they could deter- 
mine a precise cluster age with an uncertainty of only +0.3 Gyr 
for a given stellar model. However, they were unable to deter- 
mine which of the models (if any) to trust, because the difference 
in the predicted age, depending on which specific stellar model 
they adopted, was about four times greater than their measure- 
ment precision. 

In Paper I, we undertook an analysis of three detached eclips- 
ing binary systems, VI 8, V20 and V80, in NGC 6791, and deter- 
mined accurate masses and radii for the components of the first 
two of these systems. We also measured spectroscopic T e ff and 
[Fe/H] values from disentangled spectra of the binary stars and 
a few single stars. This was used, together with the cluster CMD, 
in a first demonstration of how measurements of multiple eclips- 
ing binaries in a star cluster constrain stellar models and cluster 
parameters, like age and helium content, better than previously 
possible. Here we add additional constraints from asteroseismol- 
ogy and do an extensive model comparison to determine a best- 
estimate of the helium content and cluster age while quantifying 
remaining age dependencies. As it turns out, this analysis also 
results in new insights on evolved cluster stars. 

2. Observational data 

For the present analysis we employ the observational data and 
measurements derived in paper I with slight modifications and 
additions as described in this section. 

2.1. Photometry & binary measurements 

In paper I we used the photometry of Stet son et al.l (I2003I) for 
some of the binary measurements, while the CMDs were gen- 
erated from a newer unpublished reduction of the same data 
performed by R B. Stetson. At that time we thought that differ- 
ences would be negligible and of random character. However, we 
found later that there are systematic differences between these 
data sets, due to issues w ith the photometric zero-points, as de- 
scribed in lStetsonl (I2005I) and corrected in the new reduction. 

We have therefore redone those binary measurements which 
depend on the photometry, such that our analysis is now consis- 
tently making use of the new photometric reduction. 

Only precise measurements are useful and we therefore con- 
centrate on the measurements of the eclipsing binaries V18 and 



V20, leaving out the less precise measurements of V80 from 
paper I. In Table Q] we show our measurements for these sys- 
tems from paper I with the slightly revised values of V, My and 
(m - M)y. As seen, changes to the magnitudes are small and 
the error-weighted mean apparent distance modulus remains at 
(m - M)v = 13.51 + 0.06, the value found in paper I. 

Since the following analysis will show that the reddening and 
colour-temperature relations are rather uncertain, we calculated 
the T e ff for the V20 secondary component, which does not have 
a spectroscopic T e g, under the reasonable assumption that it has 
the same distance modulus as the V20 primary component. This 
avoids the dependence on the reddening and colour-temperature 
relations which was present for this value in paper I. We have 
also chosen not to use a reddening measurement as a direct con- 
straint for the same reasons, since our only precise reddening 
measurement in paper I relied on colour-temperature relations. 

2.2. Differential reddening correction 

Several authors have noted spreads in the CMD sequence of 
NGC 679 1 that is larger than can be explained by photometric er- 
rors a lone dBedin et al.ll2008at iTwarog et al.ll201 it iPlatais et al.l 
2011). The most plausible explanation for this is differential 
reddening effects, especially since offsets from a mean fidu- 
cial are correlated with position on the sky dTwarog et al.ll2.01ll 
IPlatais et al.ll20Tlh . 

To correct for the differential reddeni ng we used a n empir- 
ical method described in detail in Mil one et al.l d2012l) . Briefly, 
we define the fiducial main sequence for the cluster and estimate 
for each star how the observed stars in its vicinity may system- 
atically lie to the red or the blue of the fiducial sequence; this 
systematic colour offset is indicative of the local differential red- 
dening. 

In Fig. [1] we compare the original (left panel) and the cor- 
rected CMD (right panel) of NGC 679 1 . It should be noted how 
many features of the CMD become narrower and more clearly 
defined after the correction has been applied, confirming that 
most of the effects of differential reddening have been removed. 
The improvement of the CMD is particularly evident for the sub- 
giant branch (SGB) and the lower red giant branch (RGB); the 
tightness along the reddening line of these stars indicate that the 
spread we see along the SGB is mainly due to differential red- 
dening. 

From star-to-star comparisons of the original and the cor- 
rected V magnitudes we can estimate star-to-star variations in 
E(B - V) and derive the reddening map in the direction of 
NGC 6791. This is shown in Fig. [2] where we have divided the 
field of view into 8 horizontal slices and 8 vertical slices and plot 
A E(B - V) as a function of the Y (upper panels) and X coordi- 
nate (right panels). We have also divided the whole field of view 
into 64x64 boxes and calculated the average A E(B - V) within 
each of them. The resultant reddening map is shown in the lower- 
left panel where each box is represented as a gray square. The 
levels of gray are indicative of the amount of differential redden- 
ing as shown in the upper-right panel. We cannot know for sure 
that this map shows only the effects of differential reddening, 
since systematic effects from the instrument and/or the reduction 
procedure may also be present in the map. However, we repeated 
the procedure using HST photometry for a smaller field within 
our FOV, and found a strong correlation with the ground-based 
results, indicating that the main effect is non-local. 

The photometry, with and without the differential reddening 
correction, is available in Tables 5 and 6 on CDS. 
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Table 1. Astrophysical data for the eclipsing binaries V18 and V20. 



V18 V20 





Primary 


Secondary 


Primary 


Secondary 


M/Mq 


0.9955 ± 0.0033 


0.9293 ± 0.0032 


1.0868 ±0.0039 


0.8276 ± 0.0022 


r/r q 


1.1011 ±0.0068 


0.9708 ± 0.0089 


1.397 ±0.013 


0.7813 ± 0.0053 


logg (cgs) 


4.3524 ± 0.0053 


4.4319 ±0.0080 


4.1840 ± 0.0078 


4.5698 ± 0.0059 


Teff (K) fl 


5600 ± 95 


5430 ± 125 


5645 ± 95 


4860 ± 125° 


V 


18.256 ± 0.020 


18.748 ± 0.020 


17.702 ± 0.020 


19.886 ± 0.020 


BC V 


-0.084 ± 0.018 


-0.123 ± 0.034 


-0.074 ± 0.018 


-0.344° 


My 


4.76 ±0.12 


5.20 ±0.16 


4.20 ± 0.09 


6.38 ±0.15° 


(m - M) v 


13.50 ±0.12 


13.55 ± 0.16 


13.51 ± 0.09 


13.51 ± 0.09° 



Notes. Measurements from paper I, slightly updated using new photometry and a new approach for V20s. Values that are different from, or were 
not included in, paper I are in boldface. 

(o) T ef! values are spectroscopic measurements, except for r cff of V20s, which is calculated assuming the same distance modulus as for V20p. 
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Fig. 1. The B,B - I CMD of NCG 6791. Left: as observed. Right: after empirical correction for differential reddening. 



2.3. Additional model constraints 

To put additional constraints on the stellar models we include in 
our model comparisons, along with the CMDs and binary mea- 
surements, the RGB mass-loss in the form of the mean mass 
difference AM = 0.09 + 0.03(random) + 0.04(systematic) Mq 
be tween the lower R GB and the red clump (RC), as measured 
bv lMiglio et al.1 d2012l) using asteroseismology. We employ only 
the mass difference because the seismic absolute masses are 
too uncertain and potentially dominated by systematic errors to 
be useful at this point. For our model comparisons we there- 
fore mark the region on the zero-age horizontal branch (ZAHB) 



with masses from that of the RGB tip down to a mass lower 
by twice the mean mass-loss (to allow for a mass-loss disper- 
sion). We also make use of our metallicity measurement from 
paper I, [Fe/H] = +0.29 + 0.03(random) + 0.08(systematic) for 
the turn-off and early SGB region, in combination with an addi- 
tional [Fe/H] measurement described in Sect. 13.4.21 
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Fig. 2. Differential reddening map of the central part of 
NGC 6791. (x,y) coordinates are in arcseconds relative to 
19:20:53.22 +37:46:18.8 (2000.0). 



3. Comparison to theoretical models 

3.1. Stellar models 

The stellar models that are used in the analysis of the binary 
stars and the cluster CMDs, were generated using a signifi- 
cantly up dated version of the University of Victoria evolution- 
ary code (IVandenBerg. Bergbusch. & Dowlerl l2()06: and refer- 
ences provided therein) with updates summarized in paper I. 
For a much more detailed description, we refer to a forthcom- 
ing study which examines the impact of varying the abundances 
of individual heavy elements (D. VandenBerg et al. 2012, sub- 
mitted). Since the settling of heavy elements is not (yet) con- 
sidered in the Victoria models, we include also a fe w models 
calculated using the DSEP code ( Fei den et al.ll201 ll and refer- 
ences therein) which includes hea vy element diffusion using the 
description of Thoul et ah ] (ff9 94) and turbulent mixing as de- 
scribed by iRichard et all (120051) . As seen in Fig. |3] the DSEP 
and Victoria codes produce nearly identical results for identical 
input physics. We can therefore be confident that it is mainly 
the effects of heavy-element diffusion we see when comparing 
Victoria and DSEP models with identical composition but in- 
cluding heavy-element diffusion in the latter. All models are cali- 
brated in order to satisfy the s olar constraints and emp loy the sur- 
face bou ndary conditions of iRrishna Swamvl ([T966). This was 
shown bv lVandenBerg et all (120081) to give a better match to the 
RGB of the solar metallicity open cluster M67 than when surface 
boundary conditions from ID model atmospheres were used. 
Solar calibration values for each model are given in Table 



3.2. General comparison procedure 

The best starting point for an analysis of binary stars is arguably 
the mass-radius diagram because it provides the most direct 
comparison between theory and observations that can be made, 



Y = 0.30, [Fe/H] = +0.30 

2 _ Solid: Victoria (t RGBr = 9.58) 
Dashed: DSEP (t RCBT = 9.63) 

Without Diffusion 



- 



Solid: Victoria (t RGBT = 9.24) 

Dashed: DSEP (t R(;Br = 9.30) 
With He Diffusion 



3.8 



3.7 



3.6 



3.5 



log T etf 

Fig. 3. Comparison between Victoria and DSEP stellar evo- 
lutionary tracks for masses and chemical abundances (as indi- 
cated) that are relevant to an investigation of NGC 6791 and the 
binary stars V18 and V20. Note that there are four tracks in the 
figure, with those without diffusion being hottest in the turn-off 
region. 

Table 2. Solar calibration values obtained for each stellar model. 



code 


element mix 


diffusion 




y," 


"MLT 


Victoria 


AGS09 fl 


Y, Tx* 


0.01323 


0.2553 


2.005 


Victoria 


GS98 fl 


Y,T^ 


0.01652 


0.2661 


2.007 


DSEP 


GS98 fl 


Y,Z 


0.01875 


0.2743 


2.139 


DSEP 


GS98" 


Y,Z,T60 C 


0.01875 


0.2743 


2.139 


DSEP 


GS98° 


Y,Z,T65 C 


0.01749 


0.2685 


2.087 



Notes. C alibrations were done u sing surface boundary conditions ac- 
cording to Krishna Swamy ( 1966). 

<n) AGS09: lAsplund et all <2009h . GS98: iGrevesse & Sauvall dl998h 
<fc) turbulence is included, details in Va ndenberg et a l . (201 2, submit- 
ted). (c> turbulence is included, details in lFeiden et al.l 1 1201 ll) . Z, and 
Yj refer to the initial value of metal and helium mass fractions. 



being independent of uncertainties in distance, reddening, and 
colour-temperature transformations . 

For the initial analysis in paper I, we scaled the relative abun- 
dances of the heavy elements recently derived for the Sun by 
Asnl und et al.l (120091) and allowed a range in [Fe/H] from +0.2 
to +0.4 with the helium mass-fraction Y as a free parameter. 
Within those parameters we were not able to find a model which 
matched the observations of the binaries and the CMD as well 
as we had hoped/expected and agreement in the mass-radius di- 
agram was not better than at the 2-3 cr level. Here we perform 
a detailed analysis, allowing also fo r the alternative abundance 
pattern of IGrevesse & Sauvall (1 19981) as well as a variation in the 
CNO abundances. 
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Interestingly, even within this expanded parameter space we 
are still unable to find a very good overall match to all observa- 
tions; models that optimise the match in the mass-radius diagram 
of the binaries are too hot to fit the T e ff s and too young to fit the 
observed CMDs. This could be a sign that our observations have 
reached a precision level which cannot be reproduced by current 
stellar models. However, it is also likely that our precision mea- 
surements have entered a regime where systematic effects are 
dominating instead of random errors. In order to avoid overin- 
terpretation we therefore accept models which match the binary 
mass and radius measurements within 3 cr, and their T^s within 
1 cr. We chose the tighter constraint on T e ff because T e ff is highly 
correlated with [Fe/H] and the distance modulus. Hence, if the 
measured r e ffS are wrong by more than 1 cr due to random er- 
rors, so are the estimated metallicity and distance modulus, and 
the model should then match a whole other constraint set. Tests 
showed that such adjustments only makes the fit worse. The T e ff 
uncertainties include a systematic contribution of 70 K (paper I) 
and we consider it unlikely that our spectroscopic temperature 
scale is different from the real one by more than the total 1 <x 
uncertainties which are ~ 100 K. 

Since we found in paper I that optimising the fit in the mass- 
radius diagram does not lead to the best overall fit, the general 
procedure adopted to find the best parameters of a given model 
was as follows: First, the age of a model with given parameters 
was adjusted to make an isochrone match the two most mas- 
sive, and therefore most age-sensitive, binary components in the 
mass-radius diagram. The corresponding isochrone in the CMDs 
would in general not match the SGB V-magnitude at the ob- 
served distance modulus, which we held fixed at the value mea- 
sured for the binaries (see Sec. 13.3. 1 I for the age uncertainty in- 
troduced by this choice); if the predicted SGB was too bright 
we increased [Fe/H] or decreased Y in the model until the corre- 
sponding isochrone matched both the most massive components 
in the mass-radius diagram and the SGB V-magnitude. We em- 
phasize that this procedure is only possible because we have de- 
termined the spectroscopic T e ffS of the binary components and 
thus the distance modulus in a way which is independent of the 
CMD and the reddening. 

E(V - I)/E(B - V) = 1.244 was chosen for all plots. This 
corresponds to an extreme lo wer limit found in th e literature by 
using the reddening law of ICardelli et al.l (Il989|). and assum- 
ing Ry = 3.1 and effective wavelengths from iMcCalll ([2004) 
for a Vega-like star. We choose this extreme limit instead of a 
more reasonable value in order to show that even allowing an 
extreme limit, we still have difficulties obtaining consistency 
among CMDs in different colour-planes. The reddening value 
for each isochrone was chosen to maximize consistency between 
colour-planes, with least weight given to the B - V colour, since 
that turned out to be the most problematic one. 



3.3. Age dependencies 

Figs. |4][8]show our measurements of the binary components in 
mass-radius and mass-r e ff diagrams and the observed V,B — V 
and V, V-/ CMDs compared to isochrones including ZAHB po- 
sitions for stellar models with different compositions and differ- 
ences in assumed physics. Here we will comment on the Figures, 
attempting to focus separately on the different factors at play and 
how the inferred age depends on them. 



3.3.1. Measurement precision 

In Fig. [4] we show a comparison of our binary measurements and 
CMDs to isochrones assuming E(B - V) - 0.14, abundances 
scaled from the solar abundance pattern of iGrevesse & Sauvall 
dl998l) to [Fe/H] = +0.35 and helium mass fraction Y = 0.30 
while neglecting heavy-element diffusion. These assumptions al- 
lowed the overall best match to the observations. Isochrones are 
shown for the best matching age and +0.5 Gyr. As illustrated, a 
very precise age of 8.3 Gyr is determined for this specific model 
with the adopted parameters. The acceptable age range seems to 
be very small when considering the mass-radius diagram as well 
as the SGB magnitude. However, the uncertainty of the binary 
T e s measurements dominate the error in the measured distance 
modulus and these two measurements are therefore not indepen- 
dent. Because of this, one would expect the best model match to 
the r e ff measurements to also be the best match to the CMD us- 
ing the measured distance modulus. In reality, we found that for 
a number of our models the best match to the SGB in the CMD 
only matched the mass-r e fj diagram at the lower 1 cr limit of 
the effective temperatures. One could therefore instead adopt the 
approach of matching the mass-r e jf diagram as the second step 
of the comparison procedure, which corresponds to moving the 
1 cr disagreement from T e ff to the distance modulus. This would 
make the inferred age younger by up to 0.4 Gyr for some models, 
which limits the precision that can be reached for these models. 
At first sight one might think that the ZAHB position would help 
solve this ambiguity but it turns out that the slight composition 
change between the results from the two fitting approaches shifts 
the ZAHB luminosity such that it matches equally well in both 
cases. Since we do not adopt this latter procedure, we will shift 
downwards our final age estimate by 0.1-0.2 Gyr, depending on 
the consistency level of the given model, and adopt +0.3 Gyr as 
our measurement precision estimate. 



3.3.2. [Fe/H] 

A comparison of the solid and the dashed isochrones in Fig. 
shows that a correlation between [Fe/H] and the helium mass 
fraction Y allows a range of acceptable model compositions 
which fit the observations equally well. With the present obser- 
vations this ambiguity in [Fe/H] and Y of the best fit is broken 
only if we trust that the models predict the colour of the red 
clump stars very accurately - something that is very question- 
able (see Sect. 13.3.6b . This means that the isochrones alone are 
not able to constrain [Fe/H], since we could produce isochrones 
with [Fe/H]= +0.20 or +0.50 that would match as well as the 
ones shown in Fig. [5] by adjusting Y along with [Fe/H]. It is 
therefore only because of our spectroscopic measurements in pa- 
per I and Sect. 13.4. 2l of this paper that we are able to constrain 
solutions to those with [Fe/H] in the range +0.29 - +0.40. 

Despite the ambiguity, we see from Fig.|5]that the age is well 
constrained, since the effect on the age from changes in [Fe/H] 
is almost exactly compensated by the corresponding change in 
Y needed to retain the match to the observations. 

The age remains well constr ained also when we assume 
abundance s scaled from those of Asplund et al.l (120091) instead 
of those o fG revesse & Sauvall (fl 998)'. as seen by comparing the 
dashed and dashed-dotted lines in Fig. [5] Any significant differ- 
ences introduced by the change in the assumed solar abundance 
pattern are compensated by the solar calibration and by a slight 
difference in helium content. Our observations therefore do not 
give preference to any particular set of solar abundances. On the 
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- age = 8.3 Gyr,[Fe/H] = + 0.35 (GS_98),Y=0.30,E(B-V) = 0. 14 
age = 7.8 Gyr,[Fe/H] = + 0.35 (GS_98),Y=0.30,E(B-V) = 0. 14 

- age = 8.8 Gyr,[Fe/H] = + 0.35 (GS_98),Y=0.30,E(B-V) = 0. 14 




Fig. 4. Measurements of the eclipsing binaries V18 and V20 (left panels) and cluster CMDs (middle and right panels) compared to 
Victoria model isochrones and ZAHB loci following the procedure described in Sect. 13.21 Un certainties shown in t he left panels are 
1 <t and 3 <x for mass and radius, and 1 cr for T e g . GS_98 denotes the solar abundance pattern of lGrevesse & Sauvalldl998l) . The thick 
part of the ZAHB c orresponds to masses with an RGB mass-loss from to twice the mean mass loss found from asteroseismology 
jlMiglio et alJl2012h . 
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- age = 8.3 Gyr,[Fe/H] = + 0.35 (GS_98),Y=0.30,E(B-V) = 0. 14 
age = 8.2 Gyr,[Fe/H] = + 0.40 (GS_98),Y=0.31,E(B-V) = 0.14 

- age = 8.5 Gyr,[Fe/H] = + 0.40 (AGS09),Y = 0.30,E(B-V) = 0. 13 




Fig. 5. As Fig. [H but sho win g isochrones of differ ent compositions. GS_98 and AGS09 denote the solar abundance patterns of 
iGrevesse & Sauvall dl998l) and lAsplund et all (120091) . respectively. 
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age = 8.2 Gyr,[Fe/H] = + 0.30 (GS_98),Y=0.30,E(B-V) = 0. 15 
age = 8.2 Gyr,[Fe/H] = + 0.35 (GS_98),Y=0.30,E(B-V) = 0. 15 




Fig. 6. As Fig.[4]but allowing for a small difference in [Fe/H] between the V18 and V20 systems, while the components of each 
system are still assumed to have identical [Fe/H]. Error boxes and ellipses are shown in the same colour as the isochrones meant to 
match them. 
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age = 8.5 Gyr,[Fe/H] = + 0.35 (GS_98),Y=0.30,E(B-V) = 0. 15 



5800 




0.9 1.0 
M/M Q 

Fig. 7. As Fig.|4]but showing a DSEP isochrone including heavy-element diffusion and turbulent mixing. In the middle panel the 
surface [Fe/H] values are shown at selected positions along the isochrone. Each star symbol marks the mean V-band magnitude 
of the stars in a spectroscopic investigation. The corresponding n umbers are th e meas u red fFe/H1 and the n u mber of stars used. I n 
order of increasing V-m agnitude, the measurem ents are those of lOriglia et al.l (12006). Car raro et al.l (120061) . iGratton et al.l d2006). 
iPeterson & GreenT(l 19981) plus the present paper, iBoesgaard et al.l (120091) . and paper I. In the right panel the change from initial t^ 
present surface [Fe/H] is shown at selected positions along the isochrone. 
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age = 8.3 Gyr,[Fe/H] = + 0.35 (GS_98),Y = 0.30,E(B-V) = 0. 14 

age = 7.5 Gyr,[Fe/H] = + 0.30 (GS_98),Y = 0.30,E(B-V) = 0. 14 
- age = 8.5 Gyr,[Fe/H] = + 0.35 (GS_98), [CNO/Fe] = -0.2,Y = 0.30,E(B-V) = 0. 14 




Fig. 8. As Fig.|4]but showing isochrones of different compositions. 
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positive side this means that the inferred age does not depend 
upon the particular choice of solar abundances. 

The comparisons so far assumed that the two eclipsing bi- 
naries have identical initial chemical composition. In Fig. [6] we 
show that if we allow for a difference of 0.05 dex in [Fe/H] be- 
tween V 1 8 and V20 then we can obtain a much better match in 
the mass-radius and mass-r e ff diagrams. A corresponding [Fe/H] 
variation in the CMDs is seen to be at the upper limit of what 
can be allowed by the main sequence colourspread, while the 
magnitudespread of the SGB remains in accordance with obser- 
vations. The observed RGB is wide enough to accommodate the 
predicted colourspread, and the observed colour variation among 
the RC stars is also better explained in this scenario. The pres- 
ence of a star-to-star variation in [Fe/H] of 0.05 dex or less could 
therefore be real. This is, in fact, the only scenario that can match 
the measurements in the mass-radius diagram to the statistically 
expected level, assuming that the measurement uncertainty es- 
timates are correct. However, since this requires an additional 
free parameter for which there is presently no conclusive evi- 
dence, and since there are many examples in the literature where 
mass and radius measurements changed at the 1 — la level as 
better observations became available, we consider this scenario 
only as a possibility until further evidence becomes available. 
The [Fe/H] measurements of the individual binary components 
are not precise enough to be used for such evidence. If accompa- 
nied by a correlated variation in Y, then [Fe/H] variations could 
even be larger as was seen in Fig. Nevertheless, as seen in 
Fig. [6] a star-to-star variation in [Fe/H] has no significant effect 
on the inferred cluster age. 

3.3.3. Heavy element diffusion 

In Fig. [7] we show an isochrone corresponding to a DSEP model 
including heavy element diffusion and turbulent mixing with a 
reference temperature log T = 6.0. As seen, the inclusion of 
heavy element diffusion in the models does not affect the inferred 
age significantly (compare to Fig. |5j. 

In the B - V plane in Fig. [7] we show the surface [Fe/H] val- 
ues at selected positions along the isochrone with initial [Fe/H] 
= +0.35. The effect of heavy element diffusion is seen to be a 
depletion of the surface [Fe/H] on the main sequence, reaching 
about 0.1 dex lower than the initial abundance at the turn-off. 
When the stars leave the main sequence, the convection zone 
deepens and brings the metals back up into the atmosphere, caus- 
ing the surface [Fe/H] to be very similar to the initial value on 
the RGB. The surface [Fe/H] becomes slightly larger than the 
initial value on the RGB because the convection zone reaches 
depths where hydrogen has been partly processed into helium. 
The model with log T = 6.0 is arguably the currently best de- 
scription of heavy element diffusion, since, with this description 
of turbulence, a model can reprodu ce abundance me asurements 
in the globular cl uster NGC 6397 (|Ko"rn et alJl2007h as well as 
other constraints dRichard et al . 2005). However, since the turn- 
off stars in NGC 6791 are cooler and have thicker convective en- 
velopes than their counterparts in the globular clusters, the tur- 
bulent mixing is much less important here. Indeed, we found 
that removing the turbulent mixing completely from the mod- 
els does not significantly change the isochrone or predicted sur- 
face [Fe/H] values. At the other extreme, when we adopt a very 
high value for the turbulence reference temperature log T = 6.5, 
the turn-off surface [Fe/H] is still depleted by 0.06 dex, in good 
agreement with the conclusion of IProffitt & Michaud ( 1 99 1J) for 
the maximum effect of turbulent mixing on diffusion. On top of 
this, we found that the surface [Fe/H] depletion is not chang- 



ing significantly with the assumed initial [Fe/H]. Therefore, we 
show on the V — I plane of Fig. [7] a rather robust prediction for 
the [Fe/H] surface abundance change along the cluster sequence, 
wh ich can be compared to observ ations . 

lAnfh onv-Twarog et aD(l2010h noted that the the lowest mea- 
sured [Fe/H] value in the literature (at that tim e) was that found 
for early SGB stars by Boesg aard et al.l (120091) while other high- 
resolution [Fe/H] measurements from stars in later evolution- 
ary stages all give higher values, as shown in Fig. Q This sug- 
gests that the signature of heavy element diffusion might al- 
ready have been observed. While large measurement uncertain- 
ties might also play a role, it is interesting that this pattern re- 
mains when we add our [Fe/H] measurement of turn-off and 
early SGB stars from paper I and that from a high-mass blue 
straggler in Sect. 13.4.21 However, more homogeneous and higher 
quality measurements are needed to show whether these indica- 
tions are real. 

Our model comparisons give us a mixed message when it 
comes to the reality of heavy-element diffusion. The effective 
temperatures on the main-sequence are higher when including 
heavy-element diffusion, partly due to a larger mixing length 
from the solar calibration and partly due to a lower surface opac- 
ity. This improves the match in the mass-r e ff diagram. But in the 
CMDs, the combined effect of higher temperatures and lower 
surface [Fe/H] shifts the isochrones to the blue. This effect is 
particularly strong in B — V, because this index is very sensi- 
tive to the surface [Fe/H] value. Since we have already chosen 
E(V - I)/E(B - V) at the extreme lower limit, we cannot choose 
a higher reddening E(B - V) in order for the isochrone to match 
the V,B - V CMD while keeping a good match to the V, V - I 
CMD. If we were able to trust the colour- temperature transfor- 
mations, this would indicate that heavy-element diffusion in real 
stars is not proceeding the way it is described in current stellar 
models. Unfortunately we cannot trust colour-temperature trans- 
formations to this level and therefore the much weaker conclu- 
sion is that either we do not properly understand heavy-element 
diffusion or the colour-temperature transformations are deficient, 
with the latter explanation being the most likely one. 

It is interesting to note that since models including heavy- 
element diffusion require a higher value of the mixing length pa- 
rameter to match the Sun, the isochrones also become hotter and 
therefore bluer on the giant branch. Because of this, the dom- 
inant effect is an overall shift in colour while the shape of the 
isochrones in the CMDs change very little, except at cooler ef- 
fective temperat ures, where uncertainties in e.g. surface bound- 
ary conditions (|VandenBerg et al. 20081). TiO molecular lines 
dGarnavich et al.lll994t ICarnev et al. 2005), and helium content 



dGirardi et al.l2007b makes the interpretation difficult. It is there- 
fore very problematic to use only the CMD to determine cluster 
parameters from an overall best fit isochrone. 

3.3.4. CNO abundances 

The initial CNO abundances in NGC 6791 are presently not well 
known, with l arge d i fferences among publi s hed measurements 
dOriglia et all l2006t ICarretta et all 120071: iPeterson & Greenl 
1998), ranging from [CNO/Fe] ~ -0.3 to solar [CO/Fe] and very 
high [N/Fe]. In addition, these measurements were performed on 
evolved stars and some of the authors find indications of extra- 
mixing processes on the RGB making it difficult to relate the 
measured values to the initial abundances. 

To investigate the effect that the assumed initial [CNO/Fe] 
has on the age of stars in NGC 6791, we computed an isochrone 
assuming [CNO/Fe] = -0.2, which we present in Fig. [8] to- 
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gether with isochrones assuming a scaled solar composition. 
Comparison of the solid and dashed lines shows that, for a 
scaled solar composition, the age decreases, temperatures in- 
crease and the SGB brightens when [Fe/H] is decreased at a fixed 
mass-radius relation. Comparing instead the solid and dashed- 
dotted lines shows a very different picture when only CNO abun- 
dances are decreased; the age increases while the SGB brightens 
without a significant temperature change on the main-sequence. 
Since a slight increase in [Fe/H] or decrease in Y is needed in or- 
der to shift the SGB magnitude in agreement with observations 
for the CNO reduced isochrone, the effect on age is an increase 
of about 0.5 Gyr. The uncertainties in predicted T e ffS caused by 
the possible effects of heavy-element diffusion and star-to-star 
variations in [Fe/H] limits the use of the measured T e ff for con- 
clusions about the reality of CNO underabundances as well as 
useful constraints to their values. 

3.3.5. The mixing length parameter 

The mixing length parameter of our models is held fixed at the 
value found for solar calibrations, and we have thus not explored 
the effects of allowing the mixing length to vary with stellar pa- 
rameters or evolutionary stage. In the follo wing we examine the 
potent ial consequences of this limitation. Trampeda ch & Steinl 
d201 lb made theoretical predictions for the variation of the mix- 
ing length as a function of T e ff and log g based on 3D convection 
simulations. To see how such predictions would affect our re- 
sults, we make use of their Fig. 6 to obtain a value of the mixing 
length for each of our four binary components and the Sun. Since 
the solar calibrated mixing l ength for our models i s high er than 
the solar value predicted by Trampeda ch & Steinl (1201 ll) we do 
not use their absolute numbers, but calculate instead their pre- 
dicted offset from the solar value for each star. Since the two 
most massive, and therefore most age-dependent binary compo- 
nents happen to have almost exactly the same predicted value 
of the mixing length as the Sun, employing such offsets would 
not affect the inferred age and helium content for a given model 
comparison. To test the effect of a varying mixing length, it is 
therefore sufficient to calculate the predicted radii of the two 
least massive binary components at the age of one of our best 
matching models using mixing lengths obtained by adding the 
predicted offsets to the solar calibrated value. The results are 
shown in Table [3] As can be seen, the effect is to lower the pre- 
dicted radii and to increase the predicted , but not by amounts 
which significantly affects the quality of the match to the obser- 
vations since changes are at the same level or smaller than the 
measurement uncertainties. 

The r e jf and thus colour of the RGB is very sensitive to the 
precise value of the mixing length. H owever, looking again at 
Fig. 6 of iTrampedach & Steinl (1201 ll) . we see that the predicted 
difference in mixing length between the most massive binary 
component and later evolutionary stages is of the order of only 
-0.02 or less. While adopting the predicted variation in mixing 
length would shift the predicted colour of the RGB, the magni- 
tude of the change is much smaller than the uncertainty of the ob- 
served RGB position. A variable mixing length would make the 
slopes of isochrones steeper on the main sequence in the CMDs, 
but the large uncertainties in the colour-temperature transforma- 
tions makes that a weak indicator as well. 

In conclusion, the relatively small consequences of 
adopting the varia t ions in mixing length predicted by 
Trampeda ch & Steinl (1201 ll) leaves us unable to distinguish be- 
tween the case of a fixed or a varying mixing length. In com- 
parison to our measurement precision, the isochrones are sim- 



Table 3. Effect of varying the mixing length parameter with 
effective temperature and surface gravity. 
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4860 


0.7813 


Measurement uncertainty(lcr) 


0.0022 
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0.0053 


Model - fixed u M lt 


0.8276 


2.005 


4837 
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Model - variable ccmlt 


0.8276 


2.265 


4889 


0.779 


change 




+0.26 
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-0.010 


V18s 
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0.9293 
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Measurement uncertainty (1 cr) 
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0.0089 


Model - fixed a MLT 


0.9293 


2.005 


5302 


0.933 


Model - variable ccmlt 


0.9293 


2.075 


5323 


0.927 


change 




+0.07 


+21 


-0.006 


Notes. For the model with fixed 


mixing 


length the 


value was 


ob- 



tained from a solar calibration. The models with variable mixing 
length has added to this va lue a difference predicted using Fig. 6 of 
Trampedach & SteI3 ( 120111) . 



ply not significantly affected by a variable mixing length, un- 
le ss variations are larger than, or different from, the predictions 
of ITrampedach & Steinl d201 lb . Measurements of more cluster 
member eclipsing binaries with lower masses could improve this 
situation in the future. 



3.3.6. Colour-temperature transformations and the reddening 

Our inability to accurately predict the relation between colours 
and effective temperatures of stars is a major uncertainty when 
comparing observations to stellar models in the CMD — in par- 
ticular at super-solar metallicities. In our model comparisons we 
have chosen colour-temperature transformations derived from 
MARCS ID atmosphere models. However, there are significant 
differences between differ ent up-to-date colour-tem perature re- 
lations, as demonstrated bv lVanden Berg et all (|2010|) . Therefore, 
the overall fairly good match between models and the observed 
CMD is likely to be more of a coincidence than a reflection of 
reality and we caution the reader that one should not trust the 
CMD colours to iden tify the most real i stic ste llar model. 

As an example, IVandenBerg et al.l (1201 Ol) noted that chang- 
ing the value of the microturbulent velocity e in the MARCS 
atmosphere models from 1 km/s to 2 km/s results in a change 
of 0.02 mag. in the B - V colour of a star with solar param- 
eters, while colours at longer wavelengths, not including the B 
filter, are not significantly affected. The direction of this shift, 
which was not reported, is to redder colours, e.g. larger values 
of B - V. In paper I we found from our spectroscopic analysis 
that e ~ 1 km/s for the turn-off stars in NGC 6791. The MARCS 
based colour-temperature relations we use, assume e = 2 km/s. If 
we demand consistency with the spectroscopic result, we would 
need to shift the predicted B - V c olours by ~ 0.02 mag to the 
blue, assuming the offset found by IVandenBerg et al.l (1201 Ol) is 
roughly independent of T e s and logg. Since Figs. [4]-[H] already 
assume an extreme lower limit for the ratio E{V - I)/E{B — V) = 
1.244, it is obvious from these figures that such a shift would 
not allow a consistent match to the CMDs in both colour planes, 
since adjusting E(B - V) to retain agreement with observations 
in B - V would make the predicted V - I colours too red. The 
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apparent inconsistency in the choice of microturbulent velocity 
exemplifies the magnitude of uncertainty pre s ent in stellar atmo- 
sphere models; Table 7 in ICasagrande et all (1201 Ol) reveals that 
the predicted B - V colour of the Sun is redder by about 0.02 
for ATLAS9 model atmospheres compared to MARCS at the 
same assumed microturbulent velocity. One might therefore sus- 
pect that using colour-temperature relations based on ATLAS9 
model atmospheres would remove the inconsistency. Even so, 
the colours remain problematic, particularly for the models in- 
cluding heavy-element diffusion. 

The adoption of empirical colour-temperature relations does 
not lead to a resolution of the proble m. We found, interest- 
ingly, t hat the empirical B - V r elation s of Sekigu chi & Fukugital 
<2000j), iRamfrez & Melende A (120051) . and iTwarog et al.1 (l2009t) 
predict almost exactly the same shape of the SGB, the turn- 
off region and the upper main-sequence on the B - V plane, 
but with colour offsets at the turn-off of -0.015, 0, and +0.025 
mag, respectively, relative to MARCS transformations. The pre- 
dicted main sequence slope is shallower than observed, becom- 
ing worse for the isochrones in cluding heavy-element di ffusion. 
The B — V transformations of ICasagrande et al.l (1201 0b match 
those of MARCS at the SGB and the turn-off, but predicts a 
much shallower main sequence than any of the other transforma- 
tions, leading to the worst match to the B - V data. At the same 
time, however, the transformations of Casagran de et al.l fcOlO) 
lead to an excellent match of the shape of the whole main se- 
quence, turn-off and SGB on the V - I plane, in particular for 
isochrones not including heavy-element diffusion. 

Because of the big uncertainty related to the B - V colours, 
our best estimate of the reddening, E(B-V), relies on V-/ rather 
than B — V. This makes the derived value of E(B - V) dependent 
on the ratio E(V-I)/E(B-V) which is not known pre cisely. For 
reasonable values around 1 .35 dHendricks et al.l2012l Tables 11- 
13) we find E(B - V) — 0.12 - 0.15 from our isochrone fits to 
V - 1 which can be made consistent also with B - V if we allow a 
systematic shift of +0.02 or more in the predicted B - V colour. 
We adopt E(B - V) = 0.14 + 0.02 as our best reddening estimate. 

The uncertainties in the colour-temperature transformations 
renders the CMDs problematic for a precise age measurement 
because they make it impossible to use the quality of the model 
fit to the CMD to select among different model assumptions. 
Luckily we can come a long way by just relying on the V-band 
magnitude of the SGB, which is the only feature of the CMDs 
which is used in our age estimates. This still relies on a bolomet- 
ric correction (BC) in the V-band, which we get from MARCS 
atmosphere models. However, since the same source is used for 
the BC of the binary components when calculating the distance 
modulus, we are only concerned about the relative error in the 
BC over a small interval in r e ff and logg, which should be very 
small. 



3.3.7. Binarity effects 

The assumption of the model comparisons we are making is that 
the binary star components evolve as single stars and are not 
significantly affected by the binarity. For binaries of lower mass 
and very short orbital periods this is often not true, with binary 
components being larger and cooler than predicted by single-star 
evolution, currently interpreted as the efficiency of convection 
being inhibited by strong magnetic fields induced by fast rota- 
tion. Although our binary stars have relatively long orbital peri- 
ods (P~15 and 20 days) one could still be concerned whether 
they are affected similarly, although to a lower extent. While 
it is presently impossible to rule this out, there are a number 



of indications which suggests that such effects are not signifi- 
cant; In our model comparisons the measured temperatures of 
the binary components are always equal to or larger than pre- 
dicted, while the effects of binarity should make them cooler. 
The stellar radii are not correlated with rotation if one assumes 
that rotation is close being to synchronized with the orbit; VI 8 
orbits slower than V20 and yet is the system where the radii 
of the components are slightly larger than predicted, while the 
opposite is true for V20. Binari ty effects should be larger at 
lower mass dClause n et al. 2009) but the isochrones match well 
both the components of V20, which have the lowest and highest 
masses measured. Finally, the two low-mass binary components 
are always matched quite well despite that our best matching 
isochrones were found using only the two most massive compo- 
nents. This suggests that binarity effects are small, if present at 
all. We thus believe that not taking into account binarity has a 
negligible effect on our age predictions in comparison to other 
uncertainties. 



3.3.8. A prolonged star-formation? 



ITwarog et al.1 (1201 ll) put forward an idea that the increased scat- 
ter in the turn-off region of CMDs of NGC 679 1 could be caused 
by a prolonged star-formation lasting up to 1 Gyr. Their argu- 
ment is based on the finding that CMDs of the inner and outer 
parts of the cluster show turn-off regions which are offset in 
colour relative to each other. The observed SGB is, however, 
much too thin to accommodate an age spread of 1 Gyr for a sin- 
gle metallicity case, see Fig. [4] The colour offset is also easily 
explained by differential reddening and goes away under such an 
assumption, cf. Sect. 12.21 Experiments with isochrones as those 
in Fig. [6] nevertheless show that under the assumption of small 
star-to-star variations in [Fe/H] it is possible to find solutions 
where the most metal-poor stars are up to 0.3 Gyr older than 
the most metal-rich while keeping the SGB thin in the CMDs. If 
accompanied by a variation from [CNO/Fe]= -0.2 for the most 
metal-poor stars to [CNO/Fe] = 0.0 for the most metal-rich then 
an age variation of as much as 0.7 Gyr would be possible with- 
out being revealed in the CMDs, as can be inferred by combin- 
ing results from Figs. [6] and [8] Although this latter scenario re- 
quires fine-tuning, we emphasize the importance of determining 
whether or not intrinsic abundance variations are present among 
the cluster stars before a prolonged star-formation can be ruled 
out. 



3.3.9. Helium abundance and inferences for AT/AZ 

As seen in Figs. |4][8] the helium mass fraction Y depends only 
weakly on the exact [Fe/H], the assumed solar abundance pat- 
tern, and [CNO/Fe], and serves in most cases to compensate the 
effect that the other assumptions have on age. For all cases within 
our spectroscopic [Fe/H] constraints we have Y = 0.30 + 0.01. 
The uncertainty in Y due to our measurement procedure, as de- 
scribed in Sec. I3.2l and l3.3. II is at the level of +0.005 or less. 

Assuming that helium enrichment follows the enrichment of 
heavy elements in a universal way and can be parametrized as 
AT/AZ, we can use our model comparisons to put constraints on 
this helium-enrichment parameter. We calculated AT/AZ for a 
subset of our best matching models both from initial solar values 
to initial values in NGC 6791, from primordial to initial solar 
values, and from primordial to initial values in NGC 679 1 . For 
the primordial value we used Y = 0.248 (e.g. iPeimbert et al.l 
l2007l:ICvburtetal.ll2008h . 



13 



K. Brogaard et al.: Eclipsing binaries in the open cluster NGC 6791 



Interestingly, when including heavy element diffusion, the 
three calculations give almost the same value, suggesting a lin- 
ear AY/AZ all the way from primordial to the high metallicity 
of NGC 6791. The exact value of AY/AZ depends on the exact 
[Fe/H] and to a lesser extent on the assumed abundance pattern, 
but the values range from ~ 1.3 to 1.5 when heavy element dif- 
fusion is accounted for, so our present best estimate for AY/AZ 
is 1.4 + 0.1. Models neglecting heavy element diffusion do not 
reproduce this result, giving a steeper slope from solar metallic- 
ity and up, rising from 1.0 (primordial to solar) to 2.0 (solar to 
NGC 6791). 

3.3.10. Best age estimate 

Taking the results of this section as a whole, we found that 
the age of NGC 6791 is very well constrained for given as- 
sumptions of physical ingredients and composition. In addi- 
tion, the predicted age is insensitive to the exact value of 
[F e/H] and whet h er one re lies on the solar abundanc e pattern 
of lAsplundeTaTI (120091) or iGrevesse & Sauvall (1 19981) . as long 
as a scaled solar composition is assumed for NGC 6791. Effects 
on age from possible binarity effects, variable mixing length, 
colour-temperature relation errors, and star-to-star variations in 
[Fe/H] and age have been shown or argued to be small or negli- 
gible. However, the combination of uncertainties related to these 
issues makes it impossible to make firm conclusions about the 
presence of heavy-element diffusion and/or non-solar [CNO/Fe] 
in NGC 6791 without further observational evidence. 

Our present best age estimate is that of a model including 
heavy element diffusion, giving more weight to the simultane- 
ous match of the effective temperatures and distance moduli of 
the binaries and the ability to produce a linear helium enrichment 
law than questionable indications from colours in the CMD. Due 
to the uncertain and conflicting results of present CNO measure- 
ments as well as their uncertain relation to the initial CNO abun- 
dances we assume for now scaled solar abundances. 

After applying an offset based on Sect. 13.3.11 our present 
best estimate for the age of NGC 6791 is thus 8.3 + 0.3 Gyr, 
where the quoted uncertainty is that originating from the fitting 
process. To this uncertainty we have to add additional contribu- 
tions, where the questionable reality of heavy element diffusion 
adds about +0.2 Gyr, while a potential prolonged star-formation 
could shift the age estimate by about +0.2 Gyr. The uncertainty 
in the CNO abundances is now the major factor in the final er- 
ror budget. Using a rough scaling of our results, our age esti- 

mate will increase b y approximately 0.7 Gyr if [CNO/Fe] 0.3 

(ICarretta et al.l 120071) are representative of the initial values for 
NGC 6791. 



3.4. New cluster insights 

3.4.1. The horizontal branch, RGB mass and mass-loss 

In Fig. [9] we show the CMD of NGC 6791 compared to one of 
our best matching isochrones and corresponding ZAHB with 
masses indicated at selected positions. As seen, the horizontal 
branch stars in the RC of NGC 6791 actually occupy a part of 
the horizontal branch which is not horizontal in CMDs with V- 
magnitude on the Y-axis. The magnitudes and the presence of 
a sloped lower envelope of the RC stars (even more evident in 
V — I CMDs in Figs.|4][8]) is evidence that they follow the upturn 
predicted by the theor etical ZAHB, excep t for a colour-shift. 

Tests with MESA dPaxton et alj201 lT) horizontal branch evo- 
lutionary tracks for masses in the range 0.8-1 2Mq predict that 



stars remain close to their ZAHB position for most of their HB 
lifetime, and then evolve quickly towards higher luminosity and 
lower r e ff corresponding to redder colours. Therefore, RC stars 
with masses above ~ O.75M0, the beginning of the ZAHB up- 
turn, should be on the ZAHB upturn or slightly redder. This sug- 
gests that our models predict too red colours, since the RC stars 
are on the blue side of the predicted ZAHB upturn. Interestingly, 
the r e ff and therefore colour at the upturn of the ZAHB, is very 
sensitive to the stellar abundances, in particular that of Oxygen, 
providing a potential indirect way of constraining the Oxygen 
abundance, if other colour uncertainties can be overcome. 

Our model ZAHB matches very well the mass-range in 
the RC which is inferred using t he RGB mass-loss of AM ~ 
0.09M© from asteroseismology (Mi a lio et all 1201 2l) . The fact 
that the absolute masses that we get on the RGB and in the 
RC are smaller than t he corresponding values from seismol- 
ogy (Basu et all 1201 it [Miglio et ai1 l2012l) is not so surprising, 
given the rather big uncertainties, and errors of systematic char- 
acter, which could be present in the seismic measurements at this 
point. 

All our acceptable model fits give a mass on the RGB at 
the V-magnitude of the RC in the range 1.15 + 0.02M Q . Our 
isochrones do not include mass-loss, so in principle the RGB 
mass could be lower, but it is unlikely that a significant mass- 
loss has occurred at the relatively low luminosities below the 
RC. 



3.4.2. Reclassification of suspected blue horizontal branch 
stars 

Stars marked with circles in Fig. [9] were selected as suspected 
horizontal branch stars in the proper motion study by Platais et 
al. (2011). As is evident, the hottest of these stars are signifi- 
cantly fainter than predicted for ZAHB stars. [Fe/H] values are 
correlated with the helium content in our best estimate models, 
and changing these two parameters shifts the ZAHB in opposite 
direction in luminosity. Therefore, we cannot find any acceptable 
model for which the ZAHB would match these stars. If we ig- 
nore this and apply an ad-hoc magnitude offset to match the blue 
stars, the ZAHB would become too faint to match the stars in 
the RC. Thus, the stars must be blue straggler stars (BSS) rather 
than blue horizontal branch (BHB) stars, unless there is some- 
thing seriously wrong with our theoretical understanding of the 
horizontal branch. 

To show that the BSS scenario is allowed by our model 
predictions, we show in Fig. [9] an extended ZAMS and main- 
sequence evolutionary tracks for models with masses of 1.4 
and 2.2Mq, the latter corresponding to twice the tum-off mass, 
which is well-defined due to our eclipsing binary measurements. 
Since the 2.2Mq track is more luminous than the stars in ques- 
tion, they could have been produced by the merger of two normal 

main-sequence stars. 

For one of the stars, known a s 2-17 (|Kinmanll 19651). m arked 
with an open circle in Fig. [9] P eterson & Greenl (119981) car- 
ried out a spectroscopic analysis. They concluded that the star 
could be either a BSS or a BHB star, but for unknown reason 
the latter option was the only one which was noted in most 
later papers w hich cited their resu lt (e.g . lOriglia et al.l G006 
ICarretta et al.l 120071: ICarraro et al.l 120061 ICarnev et al.l 12005 ' 



Gratt on et al .1 2006; Boesgaard_elalJ[2009; mentioning just a 
few). We performe d a new spectroscopic a nalysis of the spec- 
trum obtained by Peterson & Greenl (119981). For details on the 
spectrum we refer to Peterson & Greenl (11998k The analysis 
method was the same as for the binary stars in paper I, except 
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that we did not use any prior information on log g. We used 46 
Fe I and 7 Fe II lines with equivalent widths between 25 and 100 
milli-Angstrom in the wavelength range 5300-6600 Angstrom. 
Since this is a hot sta r we employed NLTE corrections from 
iRentzsch-Holml dl996h . However, our results and conclusions 
remain unchanged if we neglect the NLTE corrections, since in 
that case the inferred log g and [Fe/H] are both reduced in each 
solution by 0. 12 dex and 0.07 dex respectively. 

The analysis shown in Table |4] revealed that due to the cor- 
relation between T e ff and log g, a large range of acceptable so- 
lutions exist when no additional information is used. However, 
by combining the spectroscopic T e fi with the the V-band magni- 
tude of the star and the precise distance modulus from our binary 
star measurements, we can calculate the radius of the star, and 
therefore either the mass needed to obtain agreement with the 
spectroscopic logg, or an alternative logg, in each case. As seen 
in Table|4] the mass turns out to be higher than the turn-off mass 
in all cases, strongly indicating that this star is a blue straggler. 
Assuming that the evolution of a BSS can be approximated by 
that of a normal star, we obtain from the models in Fig.|9]a mass 
of about 1.92Mq, whereas a mass < 0.60Mq is required for the 
BHB scenario. Table [4] shows that, for the latter option there is 
no solution with a consistent value of logg. In order to reach 
consistency in even the most favorable case, an unrealistically 
large error in the distance modulus or the bolometric correction 
of ~ 0.9 mag, more than 10 times their estimated uncertainty, 
would be needed. The discrepancy is large for all BHB solutions 
and is largest for the solutions with high [Fe/H]. This strongly 
suggests that the star is not on the BHB. Furthermore, it makes it 
impossible to match a ZAHB to the star 2-17 by invoking large 
star-to-star variations in [Fe/H] and a mass-loss which inc reases 
strongly with [Fe/H], as suggested by Kalirai et al] (120071) . 

Assuming instead a mass of ~ 1 .92 Mq allows a solution 
with consistent values of log g while also being in good agree- 
ment with the other measurements of the cluster, [Fe/H] and 
photometric r e ff. The MARCS B — V colour-transformation pre- 
dicts too blue colours, but this is in accordance with the problems 
found for the isochrones in this colour plane. The best solution 
is seen to be somewhere between solutions 3 and 4 in Table |4] 
resulting in [Fe/H ] = +0. 40 + 0.08, the same value as found by 
iPeterson & Greenl (Il998l) . At the lower uncertainty-limit of this 
measurement, a significant slope of Fel with excitation potential 
appears, which is the reason we do not consider stellar models 
with [Fe/H] below +0 .3 in an y of our model comparisons. 

IPeterson & Greenl (119981) measured relative abundances 
[C/Fe], [N/Fe] and [O/Fe] at or above zero for 2-17, while these 
quantities h ave been measured to be sub-solar for g iant stars in 
NGC 6791 (ICarretta et alJl2007t lOriglia et al.ll2006l) . If this dif- 
ference is real, and not caused by large measurement errors, then 
it is another indication that star 2-17 did not follow a standard 
single-star evolution, further supporting that star 2-17 is a BSS 
and not a BHB star. 

Considering the results of this section, star 2-17 and the other 
blue s tars at the same luminosity level, suggested bv lPlatais et al.l 
d201 ll) to be BHB stars, ar e almost certainly B SSs. In fact, since 
the proper motion study bv lPlatais et al.l (1201 ll) is expected to in- 
clude all bright cluster members and yet contains only one other 
bright blue star, except for those on the EHB, it is very unlikely 
that there are any regular BHB stars present in this cluster. 

3.5. Blue straggler stars 

When comparing the positions of the blue stragglers in Fig. |9]to 
the boundaries set by the ZAMS and the TAMS, we see that they 



are skewed to wards the TAMS and slightly beyond. Following 
the results of ISills et all (120091) this is likely because the BSSs 
formed by joining the cores of two stars without much mixing 
with the material of the envelopes, and therefore the birth-line of 
the BSSs is quite far from a single star ZAMS because the mate- 
rial in the core has already been partly processed. For the same 
reason, the lifetime of a BSS is also much shorter than that of a 
normal star of the same mass. The presence of a rather high num- 
ber of BSSs (even more evident in Fig. 1 of Platais et al. (201 1) 
than in our Fig. |9]l with short main sequence lifetimes suggests 
that BSSs are continuously being produced at a high rate and 
therefore BSSs should also be present in the later evolutionary 
stages within the cluster. Indeed, some of th e stars marked with 
circles above the red clump in Fig. [9] which [Platais et all (1201 ll) 
suggested as binary HB stars could very well be evolved BSSs 
(E-BSS). Their positions and clustering ab ove the RC in th e 
CMD fit w ith this interpretatio n according to ISills et a l. (2009). 
In addition jMiglio etalj (2012) fo und a star in th e RC which, de- 
spite being a cluster member ( Stello et ai1l201 ll) . does not show 
the seismic signatur e of a RC star, and h as a seismic mass higher 
than the RC stars (ICorsaro et al.l 120121) . This star, marked with 
a square in Fig. [9] is likely to be an E-BSS in the RGB phase. 
Continued Kepler observations will provide further insights into 
these interesting stars. 

3.6. The white dwarf cooling sequence 

As an additional twist, we note that if many stars in NGC 6791 
have evolved through the BSS channel, as suggested, then the 
white dwarf (WD) cooling sequence should also harbour a sig- 
nificant fraction of E-BSSs. 

The WD cooling sequence of NGC 6791 has become known 
in the literature for its p eculiar double peak fea ture and the at- 
tempts to explain it dBedin et al.l 120051 |20 08a b). B edin et ail 
(2008b) suggest that the bright peak is due to WD binaries, 
which would produce the observed magnitude-difference be- 
tween the two peaks. However, we note that their simulated 
CMD produces the two peaks at the same colour, whereas the 
observations show the bright er peak to be b luer than the fainter 
(Fig. 2, panels b and c in Bedin et al. 2008b). 

Therefore, we speculate that the E-BSS progenitors of some 
of the WDs we observe today might play a role in the appearance 
of the bright peak of the WD luminosity function of NGC 6791. 
This could perhaps explain the colour-offset between the two 
peaks as a mass-effect; an E-BSS would be more massive and 
bluer than a normal star arriving on the WD cooling sequence 
at the same time. However, this remains speculation to encour- 
age future investigations, as it calls for extensive and quantitative 
demonstrations which are beyond the scope of this paper. The 
binary WD scenario is also not ruled out, and it is likely that a 
combination of E-BSS and binary WDs are causing the complex 
WD cooling sequence of NGC 679 1 . 



4. Summary, conclusions and outlook 

In this paper we presented an extensive comparison of stellar 
models to photometric and spectroscopic measurements of the 
eclipsing binaries V18 and V20, which are members of the old 
open cluster NGC 679 1 . 

By performing a combined stellar model comparison with 
multiple eclipsing binaries in mass-radius and mass-r e ff dia- 
grams and cluster CMDs we showed that we can constrain stellar 
models better than ever before. This allowed us to constrain the 
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Table 4. Analysis of the hot star 2-17. 



N 


[Fe/H] from 






Excitation 


A [Fe/Hf 


A(B - VV 


A(V - /)„* 


A log g c assuming 


2-17 belongs to 


M(M y 




Fel 


Fell 






slope? 








BSS(1.92Mo) rf 


BHB(0.6M o ) c 




1 


+0.54 


+0.51 


7400 


4.18 


NO 


+0.25 


-0.05 


-0.06 


+0.29 


+0.80 


3.75 


2 


+0.47 


+0.44 


7300 


4.03 


NO 


+0.18 


-0.04 


-0.04 


+0.17 


+0.67 


2.81 


3 


+0.48 


+0.47 


7200 


3.91 


NO 


+0.19 


-0.02 


-0.02 


+0.07 


+0.57 


2.25 


4 


+0.32 


+0.30 


7100 


3.68 


+l.lcr 


+0.03 


-0.02 


+0.00 


-0.13 


+0.37 


1.42 


5 


+0.26 


+0.27 


7000 


3.62 


+ 1.6(7 


-0.03 


-0.01 


+0.02 


-0.16 


+0.34 


1.31 


6 


+0.18 


+0.18 


6920 


3.58 


+2.5o- 


-0.11 


+0.00 


+0.04 


-0.18 


+0.33 


1.27 



Notes. Table shows spectroscopic solutions N = 1-6. Boldface marks for a given indicator the solutions which are preferred. Italics mark the two 
solutions between which the highest amount of consistency can be obtained. 
Measured [Fe/H] - mean cluster [Fe/H] from paper I. 

(b) Difference between the colour predicted using the spectroscopic T cff with MARCS colour-temperature relations and the observed photometric 
colour dereddened assuming E(B - V) = 0.14 and E(V - 1)1 E(B - V) = 1.244. 

(c) Difference between the spectroscopic logg and logg calculated from the spectroscopic T eff and the luminosity obtained using the observed 
V-magnitude, (m - M) v = 13.51, MARCS BC and different assumptions of the mass. 

{d) Mass for an assumed BSS star obtained from the extended ZAMS and evolutionary tracks in Fig. [9] 

w Mass for an assumed BHB star. A mass of about 0.52Mq is implied when comparing the colour of 2-17 to the model ZAHB. However, 
since in this scenario 2-17 does not lie on the ZAHB in Fig. [9] a somewhat higher mass of 0.6Mq was chosen to be sure that the mass is not 
underestimated. 

W) The mass of 2-17 derived from the spectroscopic r c(f and logg in combination with luminosity obtained from the observed V-magnitude, 
(m - M) v = 13.51, and a MARCS BC. 



helium content, and thereby obtain a more precise (and hope- 
fully accurate) age of the cluster, although in a model-dependent 
way, while demonstrating remaining age dependencies. 

Our model comparisons result in Y — 0.30+0.01, resulting in 
AF/AZ ~ 1.4, assuming that such a relation exist and applies to 
NGC 6791. It therefore seems that the cluster proto-cloud was 
rather normal despite having super-solar metallicity at an old 
age. 

NGC 6791 con tains a large number of additional detached 
eclipsing binarie s ( Rucinski et al.l Tl996: de Marchi et al. 2007; 
iMocheiskaetail |2005, Sandquist et al., in prep.). Accurate 
measurements of these, which we have shown to be possible, 
would allow even tighter model constraints by extending the 
mass range over which the observed mass-radius relation must 
be reproduced. This has the potential to strengthen the pro- 
cedure followed in this paper. Furthermore, N GC 6791 is in 
the F ield-Of-View of the NASA Kepler mission dBorucki et al.l 
1201 Oh . which will not only allow many more detached eclips- 
ing binaries to be found and measured, but also complementary 
model c onstraints from asteroseismology of the giant stars in the 
cluster dBasu et al.)l201 it iHekker et al 11201 it IPablo et al.ll201 it 
IStello et al.H201 U Icorsaro et al.l 120121; iMiglio et al.ll2012l) . All 
our acceptable model fits give a mass on the RGB at the V- 
magnitude of the RC in the range 1.15 + O.O2M0 which will 
be an important consistency check, along with the age and he- 
lium content, for future detailed asteroseismic modelling of the 
cluster giants. 

Accurate surface gravities from both eclipsing binaries and 
asteroseismology can strongly aid future spectroscopic studies 
aimed at the abundance-related issues which we have shown to 
be the dominant contributor to the remaining age uncertainty of 
NGC 6791. 

As a by-product of our analysis we were able to correctly 
classify cluster BSSs and find indications of the presence of their 
evolved counterparts. The existence of E-BSS stars in the RC 
phase can be confirmed with future Kepler data. 



Multiple detached eclipsing binaries have also been con- 
firmed in a number of o ther old open clusters (Brog aardlEoTol 
Tala mantes et alj 1201 Ol Sandquist et al., in prep.) including 
NGC6819, which is also in the Kepler Field-Of- View. Extending 
this kind of analysis to these clusters will allow more reliable 
cluster ages to be determined and more aspects of stellar models 
to be tested, when stellar models have to reproduce observations 
in clusters with a range in age and metallicity. It is our hope 
that this will also help disentangle some of the issues related to 
colour- temperature transformations and/or the model T e g scale. 

By extending studies also to the globular clusters, detached 
eclipsing binaries could ultimately provide the strongest con- 
straints on their ages and helium contents. 
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B-V 



Fig. 9. The V, B-V CMD of NGC 6791 compared to the isochrone andZAHB from Fig. |4]extended to bluer colours. Only stars with 
a non-zero proper motion probability in preliminary measurements of Kyle Cu dworth (private com m.) are shown. Circles represent 
stars which were selected as suspected HB stars in the proper motion study bv lPlatais et alj d201 ll) . The open circle, which belongs 
to this group, is the star 2-17. The square represents a clust er member that does not show the seis mic signature of a RC star and has a 
seismic mass which is larger than that of the RGB stars (Miglio "et al.ll2012llCorsaro et al .1 120121) . Also shown is an extended ZAMS 
and main-sequence evo l ution ary tracks for stars with masses of 1.4 and 2.2 Mq taking into account convective core overshooting 
as inlVandenBer g et ail (2006). The ends of the tracks are connected to show the terminal age main sequence (TAMS) in this mass 
interval. The dashed line is the 2.2 Mq track shifted down along the ZAMS to estimate the mass of 2-17. Numbers indicate masses 
at the given positions. 
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